Mixing-demixing transition and collapse of a vortex state 
in a quasi-two-dimensional boson-fermion mixture 
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We investigate the mixing-demixing transition and the collapse in a quasi-two-dimensional de- 
generate boson-fermion mixture (DBFM) with a bosonic vortex. We solve numerically a quantum- 
hydrodynamic model based on a new density functional which accurately takes into account the 
dimensional crossover. It is demonstrated that with the increase of interspecies repulsion, a mixed 
state of DBFM could turn into a demixed state. The system collapses for interspecies attraction 
above a critical value which depends on the vortex quantum number. For interspecies attraction 
just below this critical limit there is almost complete mixing of boson and fermion components. 
Such mixed and demixed states of a DBFM could be experimentally realized by varying an external 
magnetic field near a boson-fermion Feshbach resonance, which will result in a continuous variation 
of interspecies interaction. 
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I. INTRODUCTION 

A quantum degenerate Fermi gas (DFG) cannot be 
achieved by evaporative cooling due to a strong repulsive 
Pauli-blocking interaction at low energies among spin- 
polarized fermions [lj. Trapped DFG has been achieved 
only by sympathetic cooling in the presence of a second 
boson or fermion component. Recently, there have been 
successful observation 0, 0, H, |3 and associated experi- 
mental 0, i, and theoretical]!, |, M, HO, 03, El EI 
Il5l HH studies of degenerate boson-fermion mixtures by 
different experimental groups [JJ3, [H, 0| in the following 
systems: 7 Li- 6 Li Q, 23 Na- 6 Li [J and 87 Rb- 40 K %, @. 
Moreover, there have been studies of a degenerate mix- 
ture of two components of fermionic 40 K [l[ and 6 Li Q 
atoms. The collapse of the DFG in a degenerate boson- 
fermion mixture (DBFM) of 87 Rb- 40 K has been observed 
and studied by Modugno et al. 0, [l3|, [l6| , and has also 
been predicted in a degenerate fermion-fcrmion mixture 
of different-mass atoms [l7j . 

Several theoretical investigations 0, [ill, E3] of a 
trapped DBFM considered the phenomenon of mixing- 
demixing in a state of zero angular momentum when 
the boson-fermion repulsion is increased. For a weak 
boson-fermion repulsion both a Bose-Einstein condensate 
(BEC) and a DFG have maxima of probability density 
at the center of the harmonic trap. However, with the 
increase of boson-fermion repulsion, the maximum of the 
probability density of the DFG could be slowly expelled 
from the central region. With further increase of boson- 
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fermion repulsion, the DFG could be completely expelled 
from the central region which will house only the BEC. 
This phenomenon has been termed mixing-demixing in 
a DBFM. The phenomenon of demixing has drawn some 
attention lately as in a demixed state an exotic configu- 
ration of the mixture is formed, where there is practically 
no overlap between the two components and one can be 
observed and studied independent of the other. It has 
been argued @ that such a demixed state in a DBFM 
should be possible experimentally by increasing the in- 
terspecies scattering length near a Feshbach resonance 
[lcl | . More recently mixing-demixing has been studied in 
a degenerate fermion- fermion mixture [191 ]. 

On the other hand if the interspecies interaction is 
turned attractive and its strength increased, the DBFM 
collapses above a critical strength. There has been sev- 
eral experimental 0, [2(1 Hl[ and theoretical [1, EH, Ell 
studies of collapse in a DBFM of 40 K- 87 Rb mixture. As 
the interaction in a pure DFG at short distances is repul- 
sive due to Pauli blocking, there cannot be a collapse in 
it. A collapse is possible in a DBFM in the presence of 
a sufficiently strong Boson-fermion attraction which can 
overcome the Pauli repulsion among identical fermions 

It is pertinent to see how mixing-demixing manifest 
for interspecies attraction below the critical value for 
collapse. The appearance of a quantized bosonic vor- 
tex state is the genuine confirmation of superfluidity in 
a trapped BEC. In view of many experimental studies 
of such bosonic vortex states it is also of interest to see 
how the mixing-demixing phenomenon in a DBFM mod- 
ify in the presence of a bosonic vortex. The vortices 
are quantized rotational excitations [22| and can be ob- 
served in two-dimensional (2D) systems. The lowest of 
such excitations with unit angular momentum (K) per 
atom is the nonlinear extension of a well-understood lin- 
ear quantum state 23]. Vortex states in a BEC have 
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been observed experimentally (24|. Different techniques 
for creating vortex states in BEC have been suggested 
[25[, e.g., stirring the BEC with an external laser [2(|, 
forming spontaneously in evaporative cooling [27| . using 
a "phase imprin ting method" [281 ] , and rotating an axially 
symmetric trap [291 ] . Recently, the stability of the vortex 
state and the formation of persistent currents have been 
theoretically analyzed also in toroidal traps [3(J HH, . 
The generation of vortex in degenerate fermions is much 
more complex (33j and we shall not consider this possi- 
bility here. 

The purpose of this paper is to study and illustrate 
the mixing-demixing phenomenon for both attractive and 
repulsive interspecies interaction in a trapped DBFM 
vortex in a quasi-2D configuration using a quantum- 
hydrodynamic model inspired by the success of this 
model in the investigation of fermionic collapse [IB] 
and bright [U EE] and dark [2f| solitons in a DBFM. 
The conclusions of the study on bright soliton [35[ 
are in agreement with a microscopic study [37||, and 
those on collapse [l6| are in agreement with experiments 
[H, This time-dependent mean-field-hydrodynamic 

model was suggested recently [l6[ to study the collapse 
dynamics of a DBFM. 

In addition to the study of mixing-demixing in a 
DBFM in a quasi-2D configuration, we also study con- 
ditions of stability and collapse in it. Specifically, we 
study the conditions of stability when the parameters 
of a DBFM are modified, e.g., boson-boson and boson- 
fermion interactions as well as the boson and fermion 
numbers. 

There have been prior investigations of mixing- 
demixing [HI, [12] m a trapped DBFM upon an increase of 
interspecies repulsion. Also, there have been previous in- 
vestigations of stability and collapse [HI, [l5[ in a trapped 
DBFM. In contrast to these previous time-independent 
studies for stationary states, the present study relies on 
a time-dependent formulation and investigates mixing- 
demixing and stability and collapse for both attractive 
and repulsive interspecies interaction and extends to the 
case of vortex states for the first time. 

The model hydrodynamic equations in a quasi-2D form 
is derived from a Lagrangian density for the DBFM where 
the boson Lagrangian is taken in the usual mean-field 
Gross-Pitaevskii form [38|. The interaction Lagrangian 
between bosons and fermions is also taken to have the 
standard product form of boson and fermion probability 
densities. We derive a new Lagrangian for the quasi-2D 
fermions by putting them in a box of length L along x and 
y directions and in a harmonic potential well along z di- 
rection. By occupying the lowest single-particle fermion 
states we calculate the fermion probability density as a 
function of chemical potential from which we obtain the 
corresponding fermion Lagrangian density. The resultant 
hydrodynamic equations have a nonpolynomial nonlin- 
earity for the fermions, which we use in our calculation. 
In the strict 2D limit, when axial excitations are not al- 
lowed, this nonlinearity reduced to a standard cubic form. 



The paper is organized as follows. In Sec. II we present 
an account of the quantum hydrodynamic model con- 
sisting of a set of coupled partial differential equations 
involving a quasi-2D BEC and a DFG. In Sec. Ill wc 
present the numerical results on mixing-demixing and 
collapse of a DBFM in two subsections, respectively. In 
Sec. IV we present a summary and discussion. Some 
technical details are given in Appendix. 



II. BOSON-FERMION LAGRANGIAN FOR 
QUASI-2D HYDRODYNAMICS 

We consider a DBFM with Nb Bose-condensed atoms 
of mass niB and Np spin-polarized fermions of mass mp 
at zero temperature. A natural choice for a quasi-2D 
trap-geometry is a very strong confinement along the z 
axis: in this axial direction we choose a harmonic po- 
tential of frequency lu z . In the cylindric radial directions 
we take two generic external potentials for bosons and 
fermions: Vb(p) and Vf(p), where p = (x 2 + y 2 ) 1 ^ 2 is the 
cylindric radial coordinate. 

To describe this quasi-2D DBFM we use two dynam- 
ical fields: ?Pb{Pi t) and i/jf(p, t). The complex function 
ifiB{p,t) is the hydrodynamic field of the Bose gas, such 
that tib — IV'sl 2 is the 2D bosonic probability density 
and vb — id p In (i/'s/IV'sl) is the bosonic velocity. The 
complex function tl>p(p,t) is the hydrodynamic field of 
the Fermi gas, such that np = \i/jf\ 2 is the 2D fermionic 
density and vf = id p In (^f/|^f|) is the fermionic veloc- 
ity. These two complex fields are the Lagrangian vari- 
ables of the Lagrangian density 

C = L B + C F + C BF , (1) 

where Cb is the bosonic Lagrangian, Cp is the fermionic 
Lagrangian and Cbf is the Lagrangian of the boson- 
fcrmion interaction. It is important to stress that in our 
model, based on quantum hydrodynamics [38l . |39| . the 
bosonic Lagrangian C b describes very accurately all the 
dynamical properties of the dilute quasi-2D BEC [39ll40|. 
while the fermionic Lagrangian Cp can be safely used 
only for static and collective properties of the quasi-2D 
Fermi gas 39]. Note that recently quantum hydrody- 
namics has been also successfully applied to investigate 
the dimensional crossover from a 3D BEC to a ID Tonks- 
Girardeau gas [H E3] in a DBFM. 

The bosonic Lagrangian is given by [lj| HH, H(| 

= l 4(r B d^B-iPBd t rB)~7?—\V P 4>B\ 2 
2 2itlb 

h 2 l 2 

- ?n B - £b(«b) - Vbu b , (2) 

LrriBP 

where Ti 2 l 2 /(2rriBP 2 ) is the centrifugal term of the 
bosonic vortex, I is the integer quantum number of cir- 
culation and hi is the angular momentum of each atom 
in the axial (z) direction [23J. The term Eb^b) is the 
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bulk energy density of the dilute and interacting quasi- 
2D BEC under axial harmonic confinement. As shown 
in Ref. [!(| , this bulk energy density is a nonpolynomial 
function of the 2D bosonic density ns — \^b\ 2 - For small 
bosonic densities, i.e. for < tib < l/(2^/2naBBazB) 
where a sB is the 3 D Bosc-Bose scattering length and 
a Z B = V"7 \ m B^>z) is the characteristic length of axial 
harmonic confinement for bosons, the BEC is strictly 2D 
and one finds 



F - 1 2 



(3) 



where gsB = ^Ti 2 aBB / '(V^o-zbi^b) is the 2D inter- 
atomic strength [4(|. For very large densities, i.e. for 
us ^S> l/(2v / 27raBBa z s), the BEC is instead 3D and 
£b scales as n 5 J 3 (for details see Ref. [4(j, where non- 
polynomial Schrodinger equations are derived for cigar- 
shaped and disk-shaped BECs starting from the 3D 
Gross-Pitaevskii Lagrangian). Here we consider a BEC 
with a small a Z B (strong axial confinement) and a much 
smaller scattering length asB and so the BEC is strictly 
2D. 

The fermionic Lagrangian is given by [Ull 



ih 

£f = y {ip F d t ipF ~ ipFdtipp) 
— £p(np) — Vpnp , 



6lTlF 



(4) 



where £p{uf) is the bulk energy density of a non- 
interacting quasi-2D Fermi gas at zero temperature and 
under axial harmonic confinement. As shown in Ap- 
pendix, this bulk energy density is a nonpolynomial 
function of the 2D fermionic density np = IV'.fI 2 - 
For small fermionic densities, i.e. for < uf < 
1/ (2na 2 zF ), the Fermi gas is strictly 2D and one finds 
£p = ?kt> z 7r(a z i?roi?) 2 , where a Z F = Wh/JfJzmp) is the 
characteristic length of axial harmonic confinement of 
fcrmions. For large densities, i.e. for np ^> 1/ {2na 2 zF ) 1 
the Fermi gas is 3D and, as shown in Appendix, one has 
£f = (4a^ F v / 7r/3)(?ia; z )n^ 2 . Contrary to the case of 
bosons, whose dimensionality depends also on cibb (that 
is very small), for fermions it is necessary to use an ex- 
tremely small a z p and a very small number of atoms to 
have a strictly 2D configuration. This is not the case of 
real experiments and so we use the formula 



hu> z 

a lF 



for < n F a zF < ^ 



ir(n F a 2 zF ) 

J-(Ann F a 2 zF - if/ 2 + ^ for n F a 2 zF > , 

(5) 

which has been deduced in Appendix and gives the full 
2D-3D crossover of an ideal Fermi gas that is uniform in 
the cylindric radial direction and under harmonic con- 
finement in the cylindric axial direction. It is interesting 
to stress that the study of 2D-3D (and 1D-3D) cross overs 
have a long history in trapped (bosonic) atoms. There 
have been careful studies of a pair of of trapped atoms 
[43| as well as of a large number of trapped atoms [44[ ■ 



In the Appendix we consider a different type of the 2D- 
3D crossover for a large number of ideal Fermi gas atoms 
distributed over different quantum states obeying Pauli 
principle. 

In the fermionic Lagrangian of Eq. §4§ the Weiszacker 
gradient term — h 2 \\7 p iI>f\ 2 / (Giuf) takes into account the 
additional kinetic energy due to spatial variation [l4[ but 
contributes little to this problem compared to the dom- 
inating Pauli-blocking term £p{nF) [S HfjJ for a large 
number of Fermi atoms. The interaction between intra- 
species fermions in the spin-polarized state is highly sup- 
pressed due to the Pauli-blocking term and has been 
neglected in the Lagrangian Cp and will be neglected 
throughout. 

Finally, the Lagrangian of the boson-fermion interac- 
tion reads [3, l47l 



£bf 

;2. 



~9bf nB n F 



(6) 



where gsF = 2nh asF / '(niR^/27: a Z B o<zf) with asF 
the 3D Bose-Fermi scattering length and m# = 
rriBnip/ (tub + nip) the Bose-Fermi reduced mass. 

The Euler-Lagrange equations of motion of the La- 
grangian density fTj with Eqs. ©,(01), and © are given 
by 



Ih-T/IB = 



t h-^P 



n z v 2 

2m i 



2,2 



- + 



2m B p- 



:+HB{nB)+V B +gBFn F 



lpB , 

(7) 



& V? 

h 

6m p 



■ A*_f("-f) +Vp + QBFriB 



i>p, (8) 



where hb = d£B/dns — gBB n s is the bulk chemical 
potential of the strictly 2D BEC and fXp = d£p/dnp is 
the bulk chemical potential, given by Eq. (fTT)) of Ap- 
pendix, of the ideal Fermi gas in the 2D-3D crossover. 
The normalization used in Eqs. ||7J) and |8j and above is 
2TTf QC \^\ 2 pdp = N J . 

For qbf — Eq. ([7|) is the 2D Gross-Pitaevskii equa- 
tion while Eq. ([5]) is essentially a time-dependent gener- 
alization of a quasi-2D version of the time- inde pen dent 
equations of motions suggested by Capuzzi at al [14j and 
Minguzzi et al. [1 51 ] to study static and collective prop- 
erties of a confined, dilute and spin-polarized Fermi gas. 
That time-independent version was a generalization of 
the Thomas-Fermi (TF) approximation for the density 
of a Fermi gas [11 . The TF approximation for a station- 
ary Fermi gas can be obtained from Eq. ||5J) by setting the 
kinetic energy term to zero. For a large number of Fermi 
, atoms, in Eq. ([5]) the nonlinear term pp(np) is much 
larger than the kinetic energy term, hence the inclusion 
of the kinetic energy in Eq. ([8]) changes the probability 
amplitude ipp only marginally. However, inclusion of the 
kinetic energy in Eq. ((8|) has the advantage of leading to 
a probability amplitude ij)p analytic in space variable p, 
whereas the TF approximation is not analytic in p [351 ] . 

As previously discussed, the Lagrangian {T]) with |2]), 
(HI) and $Q describes a DBFM under axial harmonic con- 
finement of frequency lj z and any kind of the external 
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potentials Vb{p) and Vf(p) in the cylindric radial direc- 
tions. For our investigation of bosonic vortices in pres- 
ence of fermions we take the following radial traps 



V B (p) = V F (p) = ^ttibujIp 2 



(9) 



as in the study by Modugno et al. [1_3[ and Jezek 
et al. [45], where u>± refer to the trap frequency for 
bosons. In this way the quasi-2D mixture is achieved 
from the so-called disk-shaped configuration: the DBFM 
is confined by an anisotropic 3D harmonic potential, 
V(r) — ^rriBW±(p 2 + A 2 z 2 ), where A = uo z /uj±_ is the 
trap anisotropy. The quasi-2D configuration is appro- 
priate for studying vortices in the disk-shaped geometry 
with anisotropy parameter A> 1. 



III. NUMERICAL RESULT 

The main numerical advantage of working with Eqs. 
((7|) and ((HJ is that the calculations are much faster. In 
fact, one has to solve the two coupled differential equa- 
tions with only one space variable, p. The full 3D prob- 
lem will require an enormous computational effort. In 
our numerical simulation we consider the 40 K- 87 Rb mix- 
ture and take A = 10, ui z — 2ir x 100 Hz. We take raj as 
the mass of 87 Rb and mj? as the mass of 40 K. 

We solve numerically the coupled quantum- 
hydrodynamic equations ([7]) and ([5]) for vortex quantum 
numbers I = and I — 1 by using a imaginary-time 
propagation method based on the finite-difference 
Crank-Nicholson discretization scheme elaborated in 
Ref. [48( | . In this way we obtain the ground-state of the 
DBFM at a fixed value of the vortex quantum number 
I. We discretize the quantum- hydrodynamics equations 
([7]) and ([8]) using time step 0.0003 ms and space step 
0.02 pirn. 

The scattering length clbf is varied from positive 
(repulsive) to negative (attractive) values through zero 
(non-interacting). Note that in the experiments the scat- 
tering length clbf can be manipulated in 6 Li- 23 Na and 
40 K- 87 Rb mixtures near recently discovered Feshbach res- 
onances [181 ] by varying a background magnetic field. 



A. Mixing-demixing transition 
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FIG. 1: (Color online). DBFM with BEC having vortex quan- 
tum number I = 0. Probability densities \ipj\ 2 of bosons and 
fermions as a function of the cylindric radial coordinate p. 
Nb(= 1000) and Nf(= 120) are the numbers of bosons and 
fermions, respectively. The trap anisotropy is A = 10 The four 
panels correspond to different values of the Bose-Fermi scat- 
tering length o,bf(= 0,30 nm, 100 nm, —30 nm) and fixed 
Bose-Bose scattering length ass(= 40 nm). Note that in 
panel (d) obf is negative. 



In the first part of our numerical investigation we 
consider the mixing-demixing transition in the quasi-2D 
DBFM with vortex quantum number I = 0, 1. For a suf- 
ficiently large repulsive clbf there is demixing and for 
a large attractive asF there is mixing. If the attrac- 
tive clbf is further increased there could be collapse in 
the DBFM, which we study in detail in the next subsec- 
tion. The mixing-demixing phenomenon is quite similar 
for various boson and fermion numbers and we illustrate 
it choosing N F = 120, N B = 1000, and cibb = 40 nm. 



The results of our imaginary-time calculation with vor- 
tex quantum number I = are shown in Fig. 1, where 
we plot the probability density iV'jl 2 vs. cylindric radii 
p of the stationary boson-fermion mixture in a quasi-2D 
configuration for noninteracting, repulsive and attractive 
interspecies interaction. The probability density in Figs. 
1 and 2 is normalized to unity: 2ir J pdp = 1. In 

all cases, with o,bf > 0, because of the large nonlinear 
Pauli-blocking fermionic repulsion, the fermionic profile 
extends over a larger region of space than the bosonic 
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FIG. 2: (Color online). DBFM with BEC having vortex quan- 
tum number 1 = 1. Probability densities \tpj\ 2 of bosons and 
fermions as a function of the cylindric radial coordinate p. 
Parameters as in Fig. 1. 



one. As shown in Fig. 1, in agreement with previous 
studies in the I = state, a complete mixing-demixing 
transition is found by increasing aBF from asF = to 
aBF = 100 nm. Instead, in the case of attractive boson- 
fermion interaction (ogj? < 0) we find that the fermionic 
cloud is pulled inside the bosonic one and a complete 
overlap between the two clouds is then achieved. With 
further increase in boson-fermion interaction the system 
collapses. In Fig. 1(d), where asF — —30 nm, just below 
the critical value for collapse, we find an almost complete 
mixing between the bosonic and fermionic components. 

In the second part of the investigation we consider a 
BEC with vortex quantum number I = 1 in a DBFM 
with the same parameters of Fig. 1 . The results are dis- 
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FIG. 3: (Color online). Critical number iV£ rit of bosons vs. 
Bose- Fermi scattering length o,bf in the DBFM. The number 
Nf of fermions is fixed at 1000. Vortex quantum number I 
and Bose-Bose scattering length ass are instead varied. The 
region of collapse is indicated by arrows. 



played in Fig. 2. The non-interacting case (o,bf = 0) 
is exhibited in Fig. 2(a) with N F = 120, N B = 1000, 
clbb — 40 nm. The fermionic profile in this case is quite 
similar to that in the I = state exhibited in Fig. 1(a). 
However, the bosonic profile has developed a dip near 
origin due to the I — vortex state. In Fig. 2(b) upon 
introducing a interspecies repulsion between bosons and 
fermions a demixing has started and the fermionic wave 
function is partially pushed out from the central region 
for agp = 30 nm. This demixing has increased in Fig. 
2(c) for aBF — 100 nm. The fermionic profile in Fig. 2(c) 
for the vortex state with I = 1 is quite similar to the cor- 
responding state in Fig. 1(c) for I = 0. Finally, we find 
that an attractive boson-fermion interaction increases the 
mixing of boson and fermion components and the mixing 
is maximum for a critical value of aBF before the occur- 
rence of collapse in the DBFM. The boson and fermion 
profiles for a B F = — 30 nm just below the threshold for 
collapse is shown in Fig. 2(d). In this case the mixing is 
so perfect that the fermionic profile has developed a cen- 
tral dip near p = reminiscent of a vortex state as in the 
bosonic component. However, near p = the fermionic 
wave function tPf(p) tends to a constant value and does 
not have the vortex state behavior ipB(p) ~ p ■ This 
shows that the fermionic state is not really a vortex but 
due to mixing it tends to simulate one. 



B. Stability and Collapse 

For given values of Nb, Nf, asB and angular mo- 
mentum I, always a stable configuration is achieved for 
a repulsive (positive) aBF- However, for a sufficiently 
large attractive (negative) asF, the system collapses as 
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the overall attractive interaction between bosons and 
fermions supersedes the overall stabilizing repulsion of 
the system thus leading to instability. Also, alternatively, 
for a fixed asF, Nf, o-bb and angular momentum I, the 
system may collapse for Nb greater than a critical value 
AT£ rit . This may happen for all values of the parameters 
and a typical situation is illustrated in Fig. 3, where we 
plot A£ rit vs. clbf in different cases for Nf = 1000. One 
can have collapse in a single or both components. After 
collapse the radius of the system reduces to a very small 
value. 

For a B B < the system is stable for N < N B Tlt for 
both I = and 1 as one can see from Fig. 3, where 
we plot results for clbe — and —5 nm. The region of 
instability and collapse is indicated by arrows for N > 
Ag rit . In Fig. 3 we find that the region of stability has 
increased after the inclusion of the angular momentum 
term. Due to the stabilizing repulsive centrifugal term 
ti 2 l 2 /(2m B p 2 ) the rotating (I = 1) DBFM is more stable 
than the nonrotating (I — 0) one [49j |. 

For asB > 0, the bosons have a net repulsive en- 
ergy and the system does not collapse unless the strength 
of boson-fermion attraction |obf| is increased beyond a 
critical value. In this situation an interesting scenario ap- 
pears at fixed clef, &bb, Nf and I, as Nb is increased. 
Increasing Nb from a small value past the critical value 
Ag in \ the system collapses because the attractive in- 
teraction in the Lagrangian density ^ becomes large 
enough to overcome the stabilizing repulsions in boson- 
boson and fermion-fermion subsystems. However, when 
Nb becomes very large (Nb 3> Nf) past a second crit- 
ical value, the attractive interaction in the Lagrangian 
density ([5]) will become small compared to the overall 
repulsion of the system and a stable configuration can 
again be obtained. This is clearly illustrated in Fig. 5. 
Note that, in Eq. ©, \ip B \ 2 <x A b and \ip F | 2 oc A F 
and for a fixed total number of atoms (Nb + Nf) C-bf 
becomes large for Nb ~ Nf and small for Nb ^> Np 
or Nb <C Nf- Hence, a small number Np of fermions 
should not destabilize the stable configuration of a large 
number Nb of repulsive bosons. Similarly, a small num- 
ber Nb of bosons should not destabilize the stable con- 
figuration of a large number Np of repulsive fermions. 
This feature, that is explicitly shown in Fig. 3, should 
be quite general independent of dimensionality of space. 
It is not clear why this feature was not found in the 3D 
theoretical study of Ref. [8( . 

Next we analyze how the system moves towards col- 
lapse as the parameters of the model are changed. First 
we consider the passage to collapse as the attractive 
strength of boson-fermion interaction is increased. To 
see this we plot the root-mean-square (rms) radii prms 
of bosons and fermions vs. clbf in several cases in Fig. 
4. The rms radii remain fairly constant away from the 
region of collapse. However, as cibf approaches the value 
for collapse the rms radii decrease rapidly to a small 
value signalling the collapse. In this case both bosons 



agB = 5 nm 
N F =1000 
N B = 1000 


/ = 0, B — 
1 = 0, F ■■■■ 

/= 1, B - 
/ = 1, F 


. - 





















-20 -15 -10 -5 

a BF (nm) 



FIG. 4: (Color online). Root mean square radius prms of the 
two clouds in the DBFM as a function of the Bose-Fermi scat- 
tering length aBF- Here the number Nb of bosons (labeled 
B) is equal to the number Nf of fermions (labeled F) fixed at 
1000. The results are shown for two values of the BEC vortex 
quantum number I = 0,1. 

and fermions experience collapse simultaneously. 

We also studied the collapse for a fixed attractive asF, 
clbb and Np, while Nb is varied, to demonstrate that 
stable configuration can be attained simultaneously for 
bosons and fermions for small and large Nb, e.g. for 
N B < N F or N B > Nf- For intermediate N B there 
is collapse in either bosons or fermions or both. This is 
illustrated in Fig. 5 where we plot prms vs - log(N B ) 
for Np — 1000 and asB = 5 nm for both bosons and 
fermions for I — and 1. In both cases (I = 0, 1), as Nb 
is increased from a small value, the collapse is initiated 
near Nb = A^ rl ^ w 1000 when the radii of both bosons 
and fermions suddenly drop to a small value signalling 
a collapse in both subsystems. With further increase in 
Nb near Nb ~ 8000 the bosons pass to a stable state 
from a collapsed state and the corresponding rms radii 
increase with N B - The fermions continue in a collapsed 
state with a small radii. However, near Nb ~ 5 x 10 6 the 
fermions also come out of the collapsed state and with 
the increase of Nb the fermion radius starts to increase. 
For Nb > 10 7 , a stable configuration of the DBFM is ob- 
tained as the boson-boson and fermion-fermion repulsion 
compensates for the boson-fermion attraction. 



IV. CONCLUSION 

We have used a coupled set of quantum-hydrodynamic 
equations to study the mixing-demixing transition as well 
as stability and collapse of a trapped DBFM. The model 
equations are solved by imaginary time propagation of 
the finite-difference Crank-Nicholson algorithm. In our 
analysis the Bose-Einstein condensate is strictly 2D while 
the Fermi gas is not: for this reason we have introduced a 
new fermionic density functional which accurately takes 
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FIG. 5: (Color online). Root mean square radius prms of the 
two clouds in the DBFM as a function of log(Ns), that is 
the logarithm of the number Nb of bosons (labeled B). Here 
the number Nf of fermions (labeled F) is fixed at 1000. The 
results are show for two values of the BEC vortex quantum 
number I = 0, 1. 



justed in two cases to lead to similar strengths of the 
nonlincaritics. [It has been demonstrated, similar to the 
present Bose-Fermi mixture, that a Bose-Bose mixture 
with intraspecies repulsion and interspecies attraction 
can experience collapse 50].] But the present I = 1 re- 
sults with a bosonic vortex in a Bose-Fermi mixture have 
no analogy with in the Bose-Bose case. A slowly rotating 
Bose-Fermi mixture can have a quantized bosonic vortex 
with I — 1 with no vortex in the fermions; whereas a sim- 
ilar Bose-Bose mixture should have a I = 1 vortex in both 
the bosonic components in a stable stationary configura- 
tion. Consequently, the mean-field equations satisfied by 
the slowly rotating Bose-Fermi mixture will be distinct 
from those satisfied by a slowly rotating Bose-Bose mix- 
ture and the present results for I = 1 should be distinct 
from those for a Bose-Bose mixture. 

The present findings can be verified in experiments on 
DBFMs, specially for the vortex state, thus presenting 
yet another critical test of our quantum-hydrodynamic 
model. 



into account the dimensional crossover of fermions from 
2D to 3D. In the study of the mixing-demixing transition 
we have taken the boson-boson interaction to be repul- 
sive and the boson-fermion interaction to be both attrac- 
tive and repulsive. By considering a bosonic vortex with 
quantum number I, we have found that in both I = and 
/ = 1 cases the mixing could be almost complete up to the 
critical value of boson-fermion attraction beyond which 
the system collapses. When the boson-fermion interac- 
tion is turned repulsive, there is the mixing-demixing 
transition which is regulated by the boson-fermion re- 
pulsion. 

We have also studied the collapse in I = and I = 1 
cases. The I — 1 system is found to be more stable due 
to the centrifugal kinetic term. We have investigated in 
detail how stability is affected when the boson-boson and 
boson-fermion interaction as well as boson and fermion 
numbers are varied. 

The present analysis is based on mean-field Eqs. (O 
and |8} for the Bose-Fermi mixture, which are very sim- 
ilar in structure to those satisfied by a Bose-Bose mix- 
ture [50j . In the mean-field equations for a Bose-Bose 
mixture all nonlincarities are cubic in nature. In the 
present mean-field equations for a Bose-Fermi mixture 
apart from the intraspecies Fermi (diagonal) nonlinear- 
ity arising from the Pauli principle in Eq. |8]). given by 
Eq. (117|) . all other nonlinearities are also cubic in nature. 
In Eq. (|17p the nonlinearity is partly cubic and partly 
has a different form; whereas in the strict 2D limit this 
nonlinearity is entirely cubic in nature. Bearing such a 
similarity with the mean-field equations of the Bose-Bose 
mixture, the I = results for mixing-demixing and col- 
lapse in Bose-Fermi mixture presented here are expected 
to be similar to those of a Bose-Bose mixture provided 
that the scattering lengths and trap parameters are ad- 
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Appendix 

In this appendix we derive the zero-temperature equa- 
tions of state for an ideal Fermi gas that is uniform in 
the cylindric radial direction but under harmonic con- 
finement in the cylindric axial direction. In particular 
we obtain the chemical potential /i p and the energy den- 
sity Ep as a function of the 2D uniform radial density np 
of the Fermi gas. 

Let us consider an ideal Fermi gas in a box of length L 
along x and y axis and harmonic potential of frequency 
to z along the z axis. The total number of particles is 

N F = e(v-e ixiyi J, (10) 

' x l y ' z 

where the single particle energy reads 

h 2 (2tt) 2 /9 „ ^ , L 
= 2m F L* (t ' + $ + + 2 } ' 

Here i x , i y are integer quantum numbers and i z is a nat- 
ural quantum number. Let us approximate i x and i y by 
real numbers (see also 51]. Then 

N F = ^2 / di x di y 0(/x - e ixiyli ) . (12) 
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Setting k x = ^-i x and k y — ^j-iy, the number of particles 
can be rewritten as 

°° L 2 f 

Nf = J2nf J dk xdk y 6(ji - e kxkylz ) . (13) 



The 2D density is then 
N F 1 



L 2 4tt 2 



■ ^2 / dk x dk y 6(n - e kx k y ij • (14) 



Now we re- write the density uf in the following way 
n F = f dk x dk y 9(fi~ e k:c k v o) 



1 °° /" 

4^2 H dk x dk y 0(n- e kxkvi J . (15) 



where [ip = p — hu> z /2 is the chemical potential minus 
the ground-state harmonic energy along the z axis. 

The Eq. (fT6|) can be easily inverted and we find 



Pf = ?iw z 



2itnFa 2 Z F 



for < n^a 2 ^ < ^ 



^/47rnFa 2 F - 1 for n F a^ F > ^ . 



Equation (I17|) carries the fermionic nonpolynomial non- 
linearity in quasi-2D formulation to be used in Eqs. ([7|) 
and ([5]). In the strict 2D limit, the second term in Eq. 
(fl~5| is absent and pp — 2i:hui z nFa 2 F corresponding to 
a cubic nonlinearity. 

We can also derive the energy density £p of the Fermi 
gas from the following formula of zero-temperature ther- 
modynamics 



The first term is the density of a strictly 2D Fermi gas 
(only the lowest single-particle mode along the z axis is 
occupied) and the second term is the density which takes 
into account of all single-particle modes along the z axis, 
apart the lowest one. Setting k 2 — k 2 + k 2 the first term 
of the density np can be written as 



1 

4^2 



2 7,2 



2irkdk6(fi 



Tik 
2m 



—huj z ) = - — o 

2 ' 2na 2 



zF 



\huj z 



where a z p = yhj (myUz). The second term of the 2D 
density np is evaluated by transforming it to an integral 
and is given by 



V — 



hu z {z 27 4Tra 2 F \huj z 2 



In conclusion the 2D density np is given by 

2 



1 



np 



2ira 2 zF 



fJ-F 

frtUz 



1 / hp 

2 \ huj z 



1 



1) 



GL6) 



In this way we get 



£f = J drip /Uf(^f) 



(18) 



£ 



Tr(riFa z p) 2 for 
^n F a 2 zF -lfl 2 



< n F a 2 F < 
for 



1 

127T 



271 



n F a- zF > ± 
(19) 

The Fermi gas is strictly 2D only for < np < 
l/(2ira 2 zF ), i.e. for < fiF < ^ z . For hf > l/(2ira 2 zF ), 
i.e. for fip > fiLo z , several single-particle states of the har- 
monic oscillator along the z axis are occupied and the gas 
has the 2D-3D crossover. Finally, for uf ^> l /(2 7ra 2 P ), 
i.e. for fip ^> huj z , the Fermi gas becomes 3D |51j . 



The equations of state (JTTJ) and (fT^|) can be used to 
write down, in the local density approximation, the den- 
sity functionals of the quasi-2D Fermi gas in presence of 
an additional external potential V{p) in the cylindric ra- 
dial direction p. In this case the 2D fermionic density np 
becomes a function of the radial coordinate: uf = np{p)- 
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